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Abstract 

We study the structure and stability of the recently discussed spherically symmetric 
Brans-Dicke black-hole type solutions with an infinite horizon area and zero Hawking tem- 
perature, existing for negative values of the coupling constant lu . These solutions split into 
two classes: Bl, whose horizon is reached by an infalling particle in a finite proper time, 
and B2, for which this proper time is infinite. Class Bl metrics are shown to be extendable 
beyond the horizon only for discrete values of mass and scalar charge, depending on two 
integers m and n. In the case of even m — n the space-time is globally regular; for odd 
m the metric changes its signature at the horizon. All spherically symmetric solutions of 
the Brans-Dicke theory with uj < —3/2 are shown to be linearly stable against spherically 
symmetric perturbations. This result extends to the generic case of the Bergmann- Wagoner 
class of scalar-tensor theories of gravity with the coupling function u){4>) < —3/2. 

1. Introduction 

In the recent years there has been a renewed interest in scalar-tensor theories (STT) of gravity 
as viable alternatives to general relativity (GR), mostly in connection with their possible role 
in the early Universe. Another aspect of interest in STT is the possible existence of black holes 
(BHs) different from those well-known in GR. Thus, Campanelli and Lousto JTJ pointed out 
among the static, spherically symmetric solutions of the Brans-Dicke (BD) theory a subfamily 
possessing all BH properties, but (i) existing only for negative values of the coupling constant 
uj and (ii) with horizons of infinite area (the so-called Type B BHs These authors argued 
that large negative uj are compatible with modern observations and that such BHs may be of 
astrophysical relevance. 

In Ref. Q it was shown, in the framework of a general (Bergmann- Wagoner) class of STT, 
that nontrivial BH solutions can exist for the coupling function uj((f)) + 3/2 < 0, and that only in 
exceptional cases these BHs have a finite horizon area. In the BD theory (uj = const) such BHs 
were indicated explicitly; they have infinite horizon areas and zero Hawking temperature ("cold 
BHs"), thus confirming the conclusions of Q. These BHs in turn split into two subclasses: Bl, 
where horizons are attained by infalling particles in a finite proper time, and B2, for which this 
proper time is infinite. 

The static region of a Type B2 BH is geodesically complete since its horizon is infinitely 
remote and forms a second spatial asymptotic in a nonstatic frame of reference. For Type Bl 
BHs the global picture is more complex and is studied here in some detail. It turns out that the 
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horizon is generically singular due to violation of analyticity, despite the vanishing curvature 
invariants. Only a discrete set of Bl-solutions, parametrized by two integers m and n, admits 
a Kruskal-like extension, and, depending on their parity, four different global structures are 
distinguished. Two of them, where m — n is even, are globally regular, in two others the region 
beyond the horizon contains a spacelike or timelike singularity. 

We also discuss the stability of STT solutions under spherically symmetric perturbations. 
Under reasonable boundary conditions, it turns out that the BD solutions with u < —3/2 are 
linearly stable, and this result extends to similar solutions of the general STT provided the 
scalar field does not create new singularities in the static domain. For the case u > —3/2 the 
stability conclusion depends on the boundary condition at a naked singularity, which is hard to 
formulate unambiguously. 

Thus some vacuum STT solutions describe regular stable BH configurations with peculiar 
global structures. 



2. Black holes in scalar-tensor theories 



The Lagrangian of the general (Bergmann- Wagoner) class of STT of gravity in four dimensions 
is 



L 



-g 



<j>R + 



u(<f>) 



(1) 



where oj{4>) is an arbitrary function of the scalar field. <j) and Lm is the Lagrangian of non- 
gravitational matter. This formulation (the so-called Jordan conformal frame) is commonly 
considered to be fundamental since just in this frame the matter energy-momentum tensor Ttf 
obeys the conventional conservation law V a T° = , giving the usual equations of motion (the 
so-called atomic system of measurements) . We consider only scalar- vacuum configurations and 
put L m = 0. 

The transition to the Einstein conformal frame, g^ v = (j)~^g^ v , transforms Eq. (|lj) (up to a 
total divergence) to the form of GR with a minimally coupled scalar field ip , 



L 



'/( /'' • - 



aj3 



e = sign(w + 3/2), 



dip 



w + 3/2 



1/2 



(2) 



The field equations are 

R^v = -etp^ipv, V a V a ip = (3) 

where we have suppressed the bars marking the Einstein frame. The value e = +1 corresponds to 
normal STT, with positive energy density in the Einstein frame; the choice e = —1 is anomalous. 
The BD theory corresponds to the special case to = const, so that <fi = exp(<p/y/\u) + 3/2 1 ) . 
Let us consider a spherically symmetric field system, with the metric 

e^dfi 2 , 



dsl 



^<dt 2 



■2a 



du z 



dn 2 = d6 2 + sin 2 6d<S> 2 , 



(4) 



where E stands for the Einstein frame, u is the radial coordinate, a, /3, 7 and the field (p are 
functions of u and t. Up to the end of Sec. 3 we will be restricted to static configurations. 
The general static, spherically symmetric scalar-vacuum solution of the theory ([l]) is given 

by I! 



ds 2 



\ds% 



2bu 



J2bu 



dt 2 



du 2 



(p = Cu + (po, 



s 2 (k, u) |_s 2 (&i, u) 
C, ipo = const, 



s(k, u) 



~ l sinh ku, k > 0, 
fe = 0(j5) 



k 1 sin ku, k < 0, 



(6) 
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where J denotes the Jordan frame and u is the harmonic radial coordinate in the static space- 
time, such that a(u) = 2f3{u) + 7(11). The constants b, k and C (the scalar charge) are related 
by 



2k 2 sign k = 2b + eC" 



(7) 



The range of u is < u < u max , where u = corresponds to spatial infinity, while u max may 
be finite or infinite depending on k and the behaviour of (f>(<p). In normal STT (e = +1), by 
(0), we have only k > 0, while in anomalous STT k can have either sign. 

According to the previous studies ||, ||], all these solutions in normal STT have naked 
singularities, up to rare exceptions when the sphere u = oo is regular and admits an extension of 
the static coordinate chart. An example is a conformal scalar field in GR viewed as a special case 
of STT, leading to BHs with scalar charge j|, |6|. Even when it is the case, such configurations 
are unstable due to blowing-up of the effective gravitational coupling |7j . 

In anomalous STT (e = —1) the following cases without naked singularities can be found: 

1. k > 0. Possible event horizons have an infinite area (Type B black holes), i.e. #22 — ► 00 
as r — > 2k. In BD theory, after the coordinate transformation e 
solution takes the form 



-2ku 



1 - 2k /r = P(r) the 



dsj = P^dsl = p-t(p a dt 2 - p- a dr 2 - P^^dn 2 



Pt 



(8) 



with the constants related by (2u> + 3)£ 2 = 1 — a 2 , a = b/k. The allowed range of a and £, 
providing a nonsingular horizon at r = 2k, is 



a > 1, 



a > f > 2 



(9) 



For £ < 1 particles can arrive at the horizon in a finite proper time and may eventually (if 
geodesies can be extended) cross it, entering the BH interior (Type Bl BHs Q). When £ > 1, 
the sphere r = 2k is infinitely far and it takes an infinite proper time for a particle to reach it. 
As in the same limit 522 — ► 00, this configuration (a Type B2 BH Q) resembles a wormhole. 

2. k = 0. Just as for k > 0, in a general STT, only Type B black holes are possible 0), with 
522 ^00 as u — > co . In particular, the BD metric is 



ds 1 



-2bu 



dt 2 



C 



2bu 



du 2 



u- 



+ dn< 



s 2 (lo + 3/2) 



-2b 2 . 



(10) 



The allowed range of the integration constants is b > 0, 2b > s > —2b. This range is again 
divided into two halves: for s > we deal with a Type Bl BH, for s < with that of Type B2 
(s = is excluded since leads to GR). 

3. k < 0. In the general STT the metric @ typically describes a wormhole, with two flat 
asymptotics at u = and u = vr/|fe|, provided <j> is regular between them. In some STT the 
sphere u max = may be an event horizon, with <fi ~ l/Au 2 , Au = \u — u max | . In this case 

it has a finite area and u>(4>) + | — > —0 as it — > ii m ax- Such metrics behave near the horizon as 
the extreme Reissner-Nordstrom metric. In BD theory we have only a wormhole solution. 



For all the BH solutions mentioned, the Hawking temperature is zero. 
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3. Analytic extension and causal structure of Type Bl Brans-Dicke black 
holes 



Let us discuss possible Kruskal-like extensions of Type Bl BH metrics (|8|) and (10) of the BD 
theory. 

For (||), with a>l>£>2 — a, we introduce, as usual, the null coordinates v and w: 

v = t + x, w = t — x, x d = / P~ a dr (11) 



where x — > oo as r — > oo and x — ► —oo as r — > 2k. The asymptotic behaviour of x as r — > 2/c 
( P — > ) is x oc -P 1 " , and in a finite neighbourhood of the horizon P = one can write 

x= l (u _ u; ) = _lpl-a /( p )) (12) 

where f(P) is an analytic function of P , with /(0) = 4k /{a — 1) . Then, let us define new null 
coordinates V < and W > related to i; and u; by 

- v = (-V)-"- 1 , u; = W r - n - 1 , n = const. (13) 

The mixed coordinate patch (V,w) is defined for v < (i < — x) and covers the whole past 
horizon v = — oo. Similarly, the patch (v, W) is defined for w > (t > x) and covers the whole 
future horizon w = +oo . So these patches can be used to extend the metric through one or the 
other horizon. 

Consider the future horizon. As is easily verified, a finite value of the metric coefficient g v w 
at W = is achieved if we take n + 1 = (a — 1)/(1 — £) , which is positive for a > 1 > £. The 
metric (0) can be written in the coordinates (v, W) as follows: 

ds 2 = _ (n + 1)/ (n+2)/( n +l) . (1 _ vW n+1 )- {n+mn+1) dvdW 
Ah 2 

j-m/(n+l) . ^ _ |)r +ljm/(r.+l) f - m(i j 1 2 ( 14 ) 



(1-P) 

where m = (a — 1 + £)/(l — ■ 

The metric (|J) can be extended at fixed v from W > to W < only if the numbers n + 1 
and m are both integers (since otherwise the fractional powers of negative numbers violate the 
analyticity) . This leads to a discrete set of values of the integration constants a and £ : 

m + 1 m — n — 1 , „ 

a = , £= . 15 

m — n m — n 

where, according to the regularity conditions ([]), m > n > 0. Excluding the Schwarzschild case 
m = n + l, (£ = 0), we see that regular BD BHs correspond to integers m and n such that 

m-2>n>0. (16) 

The extension through the past horizon can be performed in the coordinates (V, w) in a 
similar way and with the same results. 

It follows that, although the curvature scalars vanish on the Killing horizon P = 0, the metric 
cannot be extended beyond it unless the constants a and £ obey the "quantization condition" 
(plf) and is generically singular. The Killing horizon, which is at a finite affine distance, is part 
of the boundary of the space-time, i.e. geodesies terminate there. A similar property was found 
in a (2+l)-dimensional model with exact power-law metric functions 
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Figure 1: Penrose diagram for a BH with even m — n 



The k = solution ( |10| ) of the BD theory also has a Killing horizon (n — ► oo) at finite 
geodesic distance if s > 0. However, this space-time does not admit a Kruskal-like extension 
and so is singular. The reason is that in this case the relation giving the tortoise-like coordinate 

x, 

lF dU= 2b^ f{u) 

(with /(oo) = 1) cannot be inverted near u = oo to obtain u as an analytic function of x. 

To study the geometry beyond the horizons of the metric (@), or fll4|), let us define the new 
radial coordinate p by 

Ik 

p = e -2ku = 1 _ m-n^ (18) 

r 

The resulting solution (|8|) , defined in the static region I (P > , p > ) , is 

d S* = P ^ dt* - 4k2 ( l™-f P— 2 dp 2 - p~™ dtf, = p—\ (19) 

By (|^), p is related to the mixed null coordinates (v,W) by 

p(v,W) = W [f{P)] 1/{n+1) [l - vW n+1 ]- 1/(n+1) ■ (20) 

This relation and a similar one giving p(V, w) show that when the future (past) horizon is 
crossed, p varies smoothly, changing its sign with W (V). For p < the metric ( |l9| ) describes 
the space-time regions beyond the horizons. Two cases must be distinguished according to the 
parity of m—n: 

1. m — n is even, i.e. P{p) is an even function. The two regions p > and p < are isometric 
(gfj,v(—p) = 9fiv{p)) f° r (l a ) n even (m even) and anti-isometric (g^u(—p) = —g^u(p)) for 
(lb) n odd (to odd). In this last case the metric tensor in region II (p < 0) has the signature 

( — h++) instead of (H ) in region I. Nevertheless, the Lorentzian nature of the space-time 

is preserved, and one can verify that all geodesies are continued smoothly from one region to the 
other (the geodesic equation depends only on the Christoffel symbols and is invariant under the 
anti-isometry — > —g^u)- The maximally extended space-time is regular, its Penrose diagram 
being the same in both cases la and lb (Fig. 1), an infinite tower of alternating regions I and 
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II. However, while in case la timelike geodesies periodically cross the horizons from one region 
to the next, in case lb these geodesies become tachyonic in region II (just as in the case of a 
Schwarzschild black hole), and generically extend to spacelike infinity (A — ► oo). 

2. m — n is odd, i.e. P(p) is an odd function; P — * — oo (r — ► 0) as p —* — oo, so that the metric 
(O) is singular on the line p = — oo which is spacelike for n odd and timelike for n even. The 
resulting Penrose diagram and geodesic motion are similar to those of Schwarzschild spacetime 
in subcase 2a (n odd, m even), and of extreme (e 2 = m 2 ) Reissner-Nordstrdm spacetime in 
subcase 2b (n even, m odd). 

4. Stability analysis 

Let us now study small (linear) spherically symmetric perturbations of the above static solutions 
(or static regions of the BHs), i.e. consider, instead of <p(u) , 

ip(u, t) = tp(u) + Scp(u, t) (21) 

and similarly for the metric functions a, (3, 7, where <p(u) , etc., are taken from the static solutions 
of Sec. 2. We are working in the Einstein conformal frame and use Eqs. (||). 

Let us choose the coordinates in the perturbed space-time so that^ 5a = 25(3 + 5j, i.e. apply 
the same coordinate condition as in the metric dss in (|5|). In this and only in this case the 
scalar equation for 5(p decouples from other perturbation equations following from (||) and reads 

e W«) S $ - 6(p" = 0. (22) 

Since the scalar field is the only dynamical degree of freedom, this equation can be used as the 
master one, while others only express the metric variables in terms of 5ip , provided the whole set 
of field equations is consistent. That it is indeed the case, can be verified directly: firstly, among 
the four different Einstein equations in (||) only three are independent and, secondly, Eq. ( p2[ ) 
may be derived from the Einstein equations. Therefore we have three independent equations for 
the three functions Sep, 5(3 and £7. 

The following stability analysis rests on Eq. (p2|). Separating the variables, 

5<p = i;(u)e iujt , (23) 

we reduce the stability problem to a boundary value problem for ip(u) . Namely, if there is a 
nontrivial solution to (^) with oj 2 < under physically reasonable boundary conditions, then 
perturbations can exponentially grow with t (instability). Otherwise it is stable in the linear 
approximation. 

Suppose — oj 2 = Q 2 , > 0. We can use two forms of the radial equation (p2[): the one 
directly following from (|23|), 

/_n2 e w^ = 0) (24) 

(as before, primes denote d/du) and the normal Liouville (Schrodinger-like) form 

d 2 y/dx 2 - [ft 2 + V(x)]y(x) = 0, V(x) = e" 4/ V - (3' 2 ). (25) 

5 The frequently used coordinate condition 5/3 = is here misleading since in our study configurations of 
utmost interest are type B black holes and wormholes, where the area function e' 3 has a minimum at some value 
of u (a throat). The equality 8(3 = is there not a coordinate condition, but a physical restriction, forcing the 
throat to be at rest. It can be shown that the condition 5/3 = creates spurious wall-like poles in the effective 
potential for perturbations, so that the latter exist separately at different sides of the throat. 
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obtained from (|24| ) by the transformation 

ip(u) = y{x) e"' 3 , x = - I e 2/3(u W (26) 



The boundary condition at spatial infinity (u = 0, x = +oo) is 5(p — > 0, or if) —> 0. For 
our metric (|5|) the effective potential ^(x) has the asymptotic form V{x) ~ 26/x 3 as n +oo, 



hence the general solution to (25) and fl24j) has the asymptotic form 



y ~ Cl e Qx + c 2 e~ nx , or ^ ~ «(ci e n/u + c 2 e- n/u ) (27) 

with arbitrary constants ci, C2- Our boundary condition leads to c\ = 0. 

For u — > Umax) where in many cases ip — ► oo, a boundary condition is not so evident. 
Refs. H 10 1 and others, dealing with minimally coupled or dilatonic scalar fields, used the 



minimal requirement providing the validity of the perturbation scheme in the Einstein frame: 

\6<p/tp\ < oo. (28) 

In STT, where Jordan-frame and Einstein-frame metrics are related by ffj^ = (l/(f>)g^ u , it 
seems reasonable to require that the perturbed conformal factor 1/0 behave no worse than the 
unperturbed one: 

\5(p/(j)\ < oo. (29) 
An explicit form (p9|) depends on the specific STT and can differ from (p8[). For example, in 



BD theory where (p oc In \<p\ , (|29|) leads to 

\5tp\ < oo. (30) 



We will refer to ( j28|) and ( |29D (or (|30|)) as the "weak" and "strong" boundary conditions respec- 
tively. If 4> and ip are regular at u — > u maiX , these conditions both give \5ip\ < oo. 

Let us discuss different cases of the STT solutions, assuming that the field 4> is regular for 
< u < n max , so that the factor in (||) does not affect the range of u. 

1. e = +1, k > 0. This is the singular solution of normal STT. As u — > oo, (3{u) ~ (b — k)u 
with b < k, so that |x| < oo and we can put x — > 0; the potential V(x) has there a negative 
double pole, V ~ — l/(4x 2 ), and (p5|) gives 

y(x) rs a/^(c3 + C4lnx). (31) 

The weak boundary condition (p8|) leads to the requirement \5(p/ip\ « |y|/iie^ w ^/(-y/xlnx) < 
oo, which is met by the general asymptotic solution ( |3l| ) and hence by its special case that joins 
the allowed solution at spatial infinity ((^7|) with c\ =0). On the other hand, in the BD theory 
the strong condition (|30|) leads to < oo as u — > oo, which is incompatible with the boundary 
condition at spatial infinity tp(u = 0) = 0, since by (24) tp" /ip > 0. 



We see that in this singular case the choice of a boundary condition is crucial for the stability 
conclusion. If we use the weak boundary condition, we conclude that the static configuration is 
unstable, in agreement with the previous work ||. This choice of the weak boundary condition 
certainly appears appropriate in the case of GR with a minimally coupled scalar field, which is 
thus unstable. On the contrary, in the BD case the strong condition seems more reasonable and 
leads to the conclusion that the BD singular solution is stable. For any other STT the situation 
must be considered separately. 
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2. e = — 1, k > 0. This case includes some singular solutions and cold BHs, exemplified by @, 
(^) for the BD theory. Now b > k, so that as u -> do, x —* — oo and V(x) ~ — l/(4x 2 ) — > 0. 
The general asymptotic solution to Eq. (^5|) has again the form (|27|). The weak condition ( p8|) 
leads, as in the previous case, to the requirement \y\/{yfx In |x|) < oo and, applied to (p7|), to 
C2 = . This means that tp must tend to zero as both u — > and u —* cxd , which is impossible 
due to ip" /ip > 0. Thus the static system is stable. Even the weak condition is sufficient to 
conclude to stability, and a stronger one is not necessary. 

3. e = — 1, k = 0. There are again singular solutions and cold BHs. Asii^oo, x — > — oo and 
the potential V(x) — > 0. The same reasoning as in item 2 leads to the same conclusion. 

4. e = — 1, /c<0. This is generically a wormhole, or in exceptional cases (see Sec. 2) a cold 
BH with a finite horizon area. In all cases one has x — > oo and y~l/|x| 3 as m -* Umax! even 
the weak boundary condition leads to \ip\ < oo, and the stability is concluded just as in items 
2 and 3. 

Thus, generically, scalar-vacuum spherically symmetric solution of anomalous STT are lin- 
early stable against spherically symmetric perturbations. Excluded are only the cases when the 
field 4> behaves in a singular way inside the coordinate range < u < u max ; such cases should 
be studied individually. 
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